Chaotic ratchet dynamics with cold atoms in a pair of pulsed optical lattices 
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We present a very simple model for realizing directed transport with cold atoms in a pair of 
periodically flashed optical lattices. The origin of this ratchet effect is explained and its robustness 
demonstrated under imperfections typical of cold atom experiments. We conclude that our model 
offers a clear-cut way to implement directed transport in an atom optical experiment. 
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I. INTRODUCTION 

The atom optics realization of the paradigmatic kicked 
rotor (KR) [1] presents the possibility to study experi- 
mentally unique quantum mechanical aspects of a fun- 
damental, classically nonlinear system. Dynamical Lo- 
calization is perhaps the most celebrated quantum phe- 
nomenon observed in the quantum KR [1, 2], but many 
other interesting features of the KR have been studied 
theoretically and experimentally [3, 4]. Very recently, 
novel applications of modified KR models have been de- 
signed which allow for a controlled, directed motion of 
particles in momentum space [5, 6]. 

The atom-optics kicked rotor (AOKR) is realized by 
subjecting cold atoms [3, 4] or a Bose condensate [7, 8, 
9, 10] to a far detuned standing wave with spatial period 
ir/k^ (fcj, being the wave number of the kicking laser) 
and pulsed with period r. The AOKR is described, in 
dimensionless units, by the Hamiltonian [11] 



H(t) 
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kcos(x)^T5(t-nT) , 



(1) 



where p is the atomic momentum in units of 2hk^ (i.e. 
of two-photon recoil momenta), x is the atomic position 
divided by 2kr,, t is time and n is an integer which counts 
the kicks. Experimentally, <5-kicks are approximated by 
pulses of width t v which are approximately rectangular 
in shape. We also define an effective Planck's constant 
ficff = T = 8E R r/h, where E R = (hk L ) 2 /2M is the recoil 
energy (acquired by an atom after emission of a photon 
with wave number k^). The dimensionless parameter 
k w VoTp/h is the kicking strength of the system (with 
Vo the height of the optical lattice creating the kicking 
potential). 

In this paper, we propose a ratchet which could be re- 
alized experimentally by adding to the standard AOKR 
dynamics defined by (1) a second kicking potential (ap- 
plied in a synchronized way with respect to the first one). 



The application of a second kicking potential to the atom 
has some analogy with the double AOKR investigated in 
[12], because in both cases a sequence of two kicks is 
applied to the atoms. In the present investigation a spa- 
tial shift of the second kick potential is also included. 
We show that this is sufficient to produce the ratchet ef- 
fect. Moreover, we consider the effects of a particle escape 
mechanism similar to evaporative cooling [13]. More pre- 
cisely, we study an open system with absorbing boundary 
conditions. If tpip) is the wave function in momentum 
space, absorbing boundary conditions are implemented 
by the prescription ip(p) = if P < ~Pc or p > p c . Such 
absorbing boundary conditions could be realized experi- 
mentally using, e.g., velocity selective Raman transitions, 
which change the internal states of the atoms, and hence 
let them escape from the states of interest [14], or by 
other state selective methods [15]. Such a scenario of 
loosing the faster atoms with momenta exceeding p c , is 
analogous to evaporative cooling of cold atoms [13]. The 
time-scale of the applied absorption mechanism should 
be of the order of the kicking period T to allow for a 
steady loss of atoms during the system's evolution. 

We point out that, as shown below, in our model the 
ratchet phenomenon is also present in the Hamiltonian 
limit without escape of particles. On the other hand, it is 
interesting to investigate the particle escape mechanism 
because it models the evaporative cooling process natural 
in cold atoms experiments. Moreover, its introduction is 
relevant in order to analyze the stability of our proposed 
ratchet mechanism after that atoms excited to higher and 
higher velocities by chaotic diffusion are eventually lost. 
Finally, particle escape may allow the unprecedented ex- 
perimental observation of a quantum phase space distri- 
bution located on an underlying classical fractal set. 

In state-of-the-art atom optics experiments, control 
over the kicking strength k (or, equivalently over the laser 
power delivered to the atoms) is achieved with a preci- 
sion of a few percent [4] . Kicking strengths in the range 
1 ... 7 correspond to standing wave amplitudes of about 



2 



80 . . . 600 Er for rubidium atoms (and assuming a rect- 
angular pulse shape with a width of 500 ns). Below we 
will be interested in the parameter region of small kick- 
ing periods T < 1, and hence it is important to note that 
time is one of the best controlled experimental param- 
eters, and kicking periods between about one hundred 
nanoseconds and a few hundred microseconds are avail- 
able, with a maximal precision of a few tens of nanosec- 
onds [4, 16, 17]. For cesium atoms, this range corre- 
sponds to dimcnsionless kicking periods T « 10~ 2 ... 20, 
and a maximal precision of 5T > 10~ 3 [17]. Atom optics 
experiments may be performed on two different atomic 
samples: laser cooled atoms and Bose-Einstein conden- 
sates. The main difference is the initial width Apo in 
momentum. For laser cooled atoms and in the best con- 
ditions, the initial width in momentum corresponds to 
a few two-photon recoils units. For Bose-Einstein con- 
densates Ap between 0.01 and 0.05 can be realized 
[7, 8, 10, 18]. Bose-Einstein condensates experience a 
nonlinear potential associated with the atom-atom inter- 
action. However, letting the condensate expand a little 
before the actual kicking evolution allows one to reduce 
the atom-atom interactions to negligible values, with only 
slight changes in Ap [7]. The present analysis focuses 
on a sample of laser cooled atoms with a large initial mo- 
mentum distribution. In fact, this condition is more fa- 
vorable for the realization of the ratchet discussed in this 
paper, because the sample explores a larger region of the 
classical phase space and therefore exploits the structure 
of phase space (a strange repeller, in the classical limit) 
induced by the evaporative cooling process. 

The paper is organized as follows: Section II analyzes 
the AOKR model and its evolution in phase space un- 
der the double kicking perturbation. Section III investi- 
gates different imperfections associated with the exper- 
imental realization. For instance a deep optical poten- 
tial is required for laser-cooled atoms, and in such condi- 
tions spontaneous emissions become a non-negligible is- 
sue. In addition, fluctuations in the laser power and other 
sources of noise are included in the analysis. The final 
Section IV concludes with an outlook discussing the role 
of nonlinearity as present in experiments using a Bose- 
Einstein condensate. 



II. MODEL AND PROPERTIES 

In this section wc introduce a kicked system that shows 
directed transport and in which the direction of the cur- 
rent can be controlled. This is done in a very simple 
way, we just have to duplicate the series of kicks in (1) in 
a convenient fashion. This simplicity is essential for an 
efficient experimental implementation with cold atoms. 

We consider a particle moving in one dimension [x £ 
(— oo, +oo)] in a periodically kicked potential. The 
Hamiltonian reads 

2 

n 2 (t) = ^- + v^{x,t), v^ = kx 



+ oo 

[S(t - riT) cos(x) + 6(t -nT - £) cos(x - <j>)\ , 

TL — — OO 

(2) 

where T is the kicking period. In fact, we propose an 
asymmetric kicking sequence. This is made out of two se- 
ries of kicks with the same spatial and temporal periods, 
2-7T and T — 2ir/ui, but shifted by a phase (j> (0 < <j> < 2ir) 
and a time £ (0 < £ < T). Due to the spatial periodic- 
ity of the kicking potential V^, the one-cycle evolution 
(Floquet) operator 

11 = g-i(r-C)p 2 /2 e -ifccos(S-0) e -i{p 2 /2 e -ifccos(£) ^ 

induced by the Hamiltonian of (2) commutes with spa- 
tial translations by multiples of 2ir. As is well known 
from Bloch theory, this implies conservation of the quasi- 
momentum f3, defined as the fractional part of the mo- 
mentum p (0 < p < 1) [19]. For a given value of 
the quasi-momentum, the wave function of the system 
is a Bloch wave, of the form e l ^ x i[)p{x), where ipp( x ) 
is a function of period 2ir. A generic wave function 
can then be written as a superposition of Bloch waves: 

V(z) = fidpeV'Mx). 

Introducing the rescaled momentum variable / = Tp, 
one can see that classical dynamics of model (2) depends 
on the scaling parameter K — kT (not on k and T sepa- 
rately). The classical limit corresponds to 7i e ff = T — > 0, 
while keeping K = fi cS k constant. 

In order to simulate the evaporative cooling process 
in the quantum model we consider the projection over a 
subspace corresponding to the quantum levels that are 
below p c (in absolute value). In practice, this is imple- 
mented at each kick: if we denote by P the projection 
operator on the interval ] — p c ,Pc[, the wave function af- 
ter n kicks is then given by 

V(p,n) = (PW) n V(p,0). (4) 

Note that quasi-momentum is still a conserved quantity. 
In the classical case, we consider lost the particles that 
reach momentum p such that \p\ > p c . 

We have checked in our numerical simulations that the 
dependence of the ratchet current on the cut-off value 
p c is weak, provided that p c 3> k. Therefore, the ratchet 
current in this regime turns out to be close to the current 
obtained in the Hamiltonian limit p c — > oo. On the other 
hand, the particle escape mechanism strongly affects the 
phase space structure, leading, in the classical limit, to 
the setting in of a strange repeller. 

In the numerical simulations reported in this paper, 
we fix K = 7, corresponding to the classically chaotic 
regime, £ = T/3, and p c h c fi = 15.2. The initial state 
is given by a uniform mixture of the momentum states 
inside the interval pU e ff £ [—1,1]. Once the quasi- 
momentum is fixed, the number of momentum states 
in this interval is oc l/7j e ff- Moreover, we average nu- 
merical data over 10 3 randomly chosen quasi-momenta. 
Classical averages are constructed from 10 7 initial con- 
ditions randomly and uniformly distributed inside the 
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region x e [0,2tt), I = pT e [-1,1]. Note that with 
these initial conditions and the above parameter values 
we are left with approximately 35% of the initial number 
of particles at time n = t/T = 10 and 10% at n = 20. 

The appearance of a strange repeller in our model in 
the classical limit is shown in the phase space portrait 
of Fig. 1 (top), obtained for = ir/2 at n — 20. The 
three panels of Fig. 1 correspond, from top to bottom, to 
the classical Poincare section and the quantum Husimi 
function at h e s ~ 0.16 and ?L e ff — 1. We can see a 
good agreement between the classical and the quantum 
phase space portraits. Quantum fluctuations smooth the 
fractal structure of the classical repeller on the scale of 
Planck's cell [20]. In the quantum case the values of 
h C H ~ 0.16 and fi. e ff — 1 considered here (and suitable 
for a realistic experimental implementation) are not suf- 
ficiently small to resolve the fractal structure at small 
scales. However, the Husimi function shows clear sim- 
ilarities with the underlying classical probability distri- 
bution. Even for h e g = T ~ 1 the major features of the 
classical repeller (i.e., width in phase space and asymme- 
try) are visible. Parameter values and evolution time are 
suitable for the experimental measurement of the quan- 
tum probability distribution located on the underlying 
classical strange repeller. This is important because the 
appearance of strange sets (attractors or repellers) is a 
distinctive feature of open chaotic systems. 

The repeller in Fig. 1 is strongly asymmetric, suggest- 
ing directed transport, that is, (p) ^ 0. This is confirmed 
by the numerical data of Fig. 2, where (p) is shown as a 
function of the time n. 

We can explain the origin of the directed current 
present in our system by following the approach devel- 
oped in [21]. We have a classical time evolution given 
by, 

x + f^(x,t) = 0, (5) 

where f$,£(x,t) — dV ( j > ^(x,t)/dx. To this equation we 
add a particle escape process consisting of cutting out 
the orbits that exceed a given value of the momentum 
p = x. We are interested in symmetry transformations 
that leave Eq. (5) invariant but change the sign of p. 
In fact, if we assume that our system is chaotic we can 
generate for each orbit its p-reversed partner, which will 
explore the whole region embedding the chaotic trajec- 
tories. This amounts to saying that, being essentially 
equivalent, both orbits (and all of them) should have 
zero average momentum. If these symmetries are ab- 
sent it is natural to conclude that a net p (i.e., different 
from zero) can be generated. Thus, breaking all possi- 
ble symmetries of this kind constitutes a good method 
to engineer ratchet systems. As the particle escape pro- 
cess introduced above is symmetrical with respect to p, 
we can neglect it in the following reasoning. It is worth 
mentioning that all the symmetry considerations devel- 
oped in this section translate almost immediately to the 
quantum case. 




FIG. 1: (Color online) Phase space pictures for cf> — 7r/2, 
at n = 20: classical Poincare sections (upper panel) and 
quantum Husimi functions at h e s ~ 0.16 (middle panel) and 
?i c ff ~ 1 (lower panel). The displayed region is given by 
/ = pT e [-20, 20] (vertical axis) and x G [0, 2ix) (horizontal 
axis). Note that, to draw the attractor, x is taken modulus 
2-7T. The color is proportional to the density: blue for zero 
and red for maximal density. 
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FIG. 2: Average rescaled momentum (7) = (p)T as a function 
of the discrete time n, for the same parameter values as in 
Fig. 1. The solid curve corresponds to the classical case, while 
the dashed curve corresponds to quantum results for Ti c g ~ 
0.16 and the dot-dashed one to fi. e tt — 1- 

There are two general ways to change the sign of p: 

(I) x — > —x + a, t — > t + 7, 

and 

(II) x — > x + a, t — > — t + 7. 

In order to leave Eq. (5) unchanged we need that 
fcf,^(x,t) = —f < f,^(—x + a,t + 'y) holds for (I), since 
x — » —x under this transformation. If we apply twice 
transformation (I) we obtain f<f>,£(x + a/2, t) = f^.^x + 
a/2,t + 2-f). Since f^,^(x,t) is assumed to be bounded 
and periodic with zero mean, both in x and t, 7 can only 
be an integer multiple of T/2 (including the 7 = case). 
In turn, there are no restrictions on a. On the other 
hand, for (//) we need (x, t) — +f ( j ) ^(x + a, — t + 7) 
(with a plus sign since now x keeps its original sign) . By 
applying twice transformation (II) we obtain /^(rr, t + 
7/2) = f ( f,^(x+2a, t+7/2). Following the same reasoning 
as before, a is fixed to integer multiples of ir (including 
a = 0) while there are no restrictions on 7. Note that 
(I) and (II) are the only two symmetries that should be 
broken in order to find directed transport. Our choice of 
the potential (2) guarantees the possibility to break both 
of them. 

In fact, we have that f^(x,t) = fcE-^[ — 5(t — 
nT) sin(x) — 5(t — nT — £) sin(x — </>)], and in the case 
of symmetry (I) we require that f^^(x,t) — — x + 

a,t + 7). We can take 7 = without loss of gener- 
ality since we only have a sum of delta functions in 
t, i.e. the sign change of / induced by symmetry (I) 
can only come from the first part of the transformation 
(x — > —x + a). Therefore, we arrive at the conditions 
sin(— x+a) = — sin(x) and sin(— x+a — <f>) = —sm(x—(j)). 
These two conditions lead to a = I2n and a — l'2n + 2<f>, 
with I and I' integers, and cannot be fulfilled together, 
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FIG. 3: Average momentum (p) as a function of n, for 
cj> — tt/2 (positive values), <f> — (zero values), and (f) = —n/2 
(negative values) Both the classical (solid curves) and the 
quantum case (dot-dashed curves, h c g ~ 1) are shown. Note 
that at 4> = quantum and classical curves are almost super- 
imposed. 

except for <\> = or <p — ft- Therefore, symmetry (I) is 
broken when <p ^ 0, tt. 

In the case of symmetry (II), if we take a an odd 
multiple of n then the sign of changes. Then, we 
are only left with a being an even multiple of n, i.e., we 
can take a = without loss of generality. Moreover, we 
notice that if (f> = and a = both kicks become the 
same in x and therefore symmetry (//) holds for any £, 
taking 7 = £. On the other hand, considering 0^0 
we arrive at the conditions J2n=-cv + 7 ~ n T) = 
E^oo *(* - nT) and £^-00 5(-t + 7 - nT £) = 

ErtT-oo^ - 0. which im P!y 7 - IT and 7 - I'T + 
2£, with Z and V integers. We conclude that, if <f> ^ 0, 
symmetry (II) is broken when £ ^ 0,T/2. 

In summary, both symmetries (I) and (II) are broken 
for ^ 0,7r and £ 7^ 0,T/2. Hence two series of kicks 
are sufficient to observe the ratchet effect, provided that 
these kicks are shifted both in space and in time, the 
shift in space being different from half wave length and 
the shift in time being different from half period. 

It is interesting to remark that current reversal can be 
engineered in a very simple way, by taking <f> = —<j) in- 
stead of 4> m (2). Indeed, Eq. (5) is left unchanged when 
x — > —x, t —* t, and — > <j> — —cf), while this transforma- 
tion changes the sign of p. We can see current inversion 
in Fig. 3, both in the classical and in the quantum case, 
when cj) = n/2 —> <f> = -tt/2. Note that (p) = at cj) = 0, 
in agreement with the above symmetry considerations. 

III. STABILITY OF THE RATCHET EFFECT 
UNDER IMPERFECTIONS 

The purpose of this section is to study the robustness of 
the ratchet effect introduced in this paper in the presence 
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FIG. 4: Average momentum (p) as a function of the discrete 
time n, for the same parameters as in Fig. 2, at h c g ~ 1. At 
each kick, the quasi-momentum can jump to any other possi- 
ble value with probabilities (dot-dashed curve), 0.2 (dashed 
curve), and 0.5 (long dashed curve). The solid curve corre- 
sponds to the classical case. 

of typical sources of noise in cold-atom experiments. For 
the large kicking strengths needed to guarantee clear sig- 
natures of a chaotic repeller, spontaneous emission dur- 
ing the flashing of the optical lattice cannot be ruled out 
[3]. Spontaneous emission can be effectively modeled by 
random jumps in quasi-momentum [19]. We test the in- 
fluence of such random changes in quasi-momentum on 
the results presented in the previous section. That is to 
say, we repeat the previous calculations but letting at 
any kick the quasi-momentum randomly change with a 
probability of 0, 0.2 and 0.5 (see Fig. 4). In practice, it 
may jump to any possible value in the Brillouin zone with 
those probabilities. As can be seen, this additional ran- 
domness even helps to reduce fluctuations, and when the 
jump probability is different from zero there is a better 
convergence towards the classical result. 

We now investigate how different kind of errors affect 
the value of the ratchet current. More precisely, we com- 
pute the average current (p) av , obtained after averaging 
(p) in the time interval 10 < n < 20, as a function of the 
noise strengths associated to different noise sources. 

First of all, we consider the effects of fluctuations in 
the kicking strength. This is simulated by memory less 
random errors of size SK in the value of K: the kicking 
strength K n at time n is given by K n = K + (5K) n where 
the noise value (8K) n is randomly drawn from a uniform 
distribution in the interval [-8K, SK}. It can be seen in 
Fig. 5 that the ratchet effect is stable up to approximately 
SK ~ 2, corresponding to a relative amplitude noise of 
SK/K ~ 0.3. 

Since the ratchet mechanism described in the pre- 
vious section works the better the smaller we choose 
h c H = T, we consider possible fluctuations in the kicking 
period [16] arising from the problem of controlling strong 
but narrow pulses in time with a high repetition rate. 



FIG. 5: Average current (p) a v as a function of the noise SK 
in the kick strength K, for parameter values as in Fig. 2, at 
heH ^ 1. 




FIG. 6: Average current (p)av as a function of the relative 
error ST /T in the kicking period T, for parameter values as 
in Fig. 5. 



We model these imperfections as random and memory- 
less fluctuations in the period between consecutive kicks. 
This takes into account the fact that the timing of the 
kicks can suffer from uncontrollable variations. As we 
can see from Fig. 6, stability is quite satisfactory when 
ST IT < 0.5, where ST is the size of the fluctuations and 
T ~ 1. 

Finally, we consider the effect of an imprecision in 
the p c selection. This is modeled by random memory- 
less variations of the cut-off value (p c )n used at time n: 
(p c )n = Pc+ {Sp c )n, with {8p c )n e [-5p c ,Sp c \. Again 
the ratchet effect proves to be robust, as can be deduced 
from Fig. 7. The results of this figure arc in agreement 
with the previous observation that the dependence of the 
ratchet current (p) on the cut-off value p c is weak (under 
the condition p c ^$> k). 
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FIG. 



7: Average current (p)av as a function of the relative 
error Sp c /p c in the cut-off momentum p c , for parameter values 
as in Fig. 5. A magnification of the figure for small values of 
Spc/pc is shown in the inset. 



IV. CONCLUSIONS 

Considering a realistic experimental scenario, we 
showed that a ratchet effect - induced by a combination 
of a two-kick sequence as applied to an open system - 
is observable in an atom-optics kicked rotor experiment. 
We also checked the robustness of the ratcheted atomic 
evolution under reasonable noise conditions. 

An interesting perspective would be to study the 
ratchet dynamics in a kicked Bose-Einstein condensate. 
Strong kicks may however, lead to thermal excitations 
out of equilibrium and destroy the condensate, rendering 
the description by the usually applied Gross-Pitaevskii 
equation meaningless [22]. We have verified that the 
ratchet evolution is preserved in the presence of typical 
experimental nonlincarities. However, a full treatment 
of a strongly kicked Bose-Einstein condensate remains a 
challenge for future work. 
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